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AN INLET ANALYSIS FOR 
THE NASA HYPERSONIC RESEARCH ENGINE 
AEROTHERMODYNAMIC INTEGRATION MODEL 


By Earl H. Andrews, Jr,, James W. Russell,* 

Ernest A. Mackley, and Ann L. Simmonds 
Langley Research Center 

SUMMARY 

An inlet theoretical analysis has been conducted in support of the aerothermodynamic 
research program for the NASA Hypersonic Research Engine (HRE) Project which is 
scheduled to culminate in an experimental investigation of the Aerothermodynamic 
Integration Model (AIM)^ in the Lewis Hypersonic tunnel facility (HTF) at the Plum Brook 
Station at nominal Mach numbers of 5, 6, and 7. The AIM experimental investigation is 
for the purpose of determining the net aerodynamic and thermodynamic effects of the 
three full-scale components (inlet, combustor, and nozzle) when integrated. It is desir- 
able to evaluate the AIM combustor performance; however, the combustor entrance condi- 
tions (inlet throat conditions) will not be experimentally obtained from AIM test measure- 
ments. Therefore, it was considered necessary to conduct a theoretical analysis of the 
AIM inlet. A method -of -characteristics computer program was employed to perform the 
inlet computations (the program can compute for two-dimensional or axisymmetric, real- 
gas or ideal -gas, supersonic external and internal flows, with or without boundary layer 
at constant wall temperatures) . 

The computer program was used for generating a number of cases for various real- 
gas, viscid-inviscid flow conditions. Ranges of mass -flow ratios and additive drag coef- 
ficients for the scheduled AIM test conditions were calculated, and "curve fits" for the 
mass -flow ratios as a function of cowl position and Mach number were obtained, Mass- 
weighted average inlet performances — total-pressure recovery, kinetic energy efficiency, 
and throat Mach number - were also calculated. "Shock -on-lip" total-pressure-recovery 
results supplemented previously documented sparse HRE full-scale (flight condition) ana- 
lytical results to the extent that a nearly straight-line decreasing trend was indicated 
over the free -stream Mach number range from 4 to 8. 


*LTV Aerospace Corporation, Hampton, Va. 

^Inlet analysis performed in 1972; AIM tests completed May 1974. 


INTRODUCTION 


As part of the NASA Hypersonic Research Engine (HRE) Project, tests of a complete 
engine (inlet-combustor -nozzle) designated the Aerothermodynamic Integration Model 
(AIM) are planned. The AIM is a full-scale, water-cooled, boilerplate engine which burns 
hydrogen fuel. Component tests have been performed on (1) a 2/3-scale inlet model 
in unheated-air wind tunnels (ref. 1); (2) a two-dimensional combustor model with direct- 
connect hot vitiated airflow (ref. 2); and (3) a nozzle model with direct -connect cold and 
hot airflow (ref. 3). Another two-dimensional combustor model, which closely repre- 
sented a segment of the HRE -AIM annulus combustor, was tested in 1971 (ref. 4). The 
experimental investigation of the AIM is scheduled to be conducted in the Lewis hypersonic 
tunnel facility (HTF) at the Plum Brook Station at nominal Mach numbers of 5, 6, and 7 
for the purpose of determining the net aerodynamic and the thermodynamic effects of the 
three components when integrated as a complete full-scale engine. 

Inlet throat (combustor entrance) conditions must be known in order to properly 
evaluate the combustor performance during the AIM experimental investigation. However, 
direct measurements for determining inlet throat conditions and performance will not be 
obtained during the AIM tests; surface pressure measurements and calculated skin- 
friction forces will be used to determine the inlet momentum losses and thus the per- 
formance. Therefore, it was considered necessary to theoretically analyze the AIM inlet 
in order to obtain throat conditions needed for evaluations of the AIM combustor perform- 
ance and inlet performance values for comparison with those to be obtained during the 
AIM tests, and to thoroughly map the inlet mass ratios for the AIM test conditions. The 
need for theoretical inlet mass -flow data became increasingly important when the engine 
airflow metering tests originally planned had to be abandoned. Some theoretical analyses 
at Mach 4,6, and 8 flight conditions were performed for the HRE full-scale inlet and 
reported in reference 1. The full-scale inlet results of the reference 1 analysis are very 
sparse and do not well define the range of conditions for the AIM experimental investiga- 
tion. Therefore, the purpose of the present analysis is to supplement the results of ref- 
erence 1 so as to have a more thorough inlet analysis for the conditions of the AIM 
experimental investigation. Such an analysis was performed by using an updated version 
of the method-of -characteristics computer program described in reference 5. The pro- 
gram was used for defining numerous inlet real -gas, viscid -inviscid flow -field conditions 
at Mach 5.0, 5.15, 6.0, 6.15, 7.0, 7.15, and 7.25 for various cowl lip positions relative to 
the inlet spike tip. Analysis results are presented in terms of schedules of mass -flow 
ratios, aerodynamic contraction ratios, and additive drag, internal-flow shock patterns, 
and performance parameters including mass -weighted average total -pressure recoveries, 
throat Mach numbers, and kinetic energy efficiencies. 
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SYMBOLS 


Values are given in both SI and U.S. Customary Units. The measurements and 
calculations were made in U.S. Customary Units. 

A area, m^ (ft2) 

A 

—22 contraction ratio 

Ath 

aerodynamic contraction ratio 

AthWooy 

Cd ^ additive drag coefficient 

h enthalpy, J/kg (Btu/lb) 

h' throat height (inclined 95.645° to the AIM center line) 

m/moo ratio of captured mass flow to that free- stream mass flow that 
passes through an area equal to the projected cowl area 
(Rc^ = (22.934 cm)2 = (9.029 in.)2 = (0.752 ft)^) 

M Mach number 

p pressure, atm (psia) 

q dynamic pressure, atm (psf) 

R radius nondimensionalized by (22.86 cm (9.0 in.)) 

T temperature, K (°R) 

V velocity, m/sec (ft/sec) 

X distance from spike vertex nondimensionalized by R^^^ (see fig. 1(a)) 

y/h’ ratio of distance from centerbody surface to throat height 

?7ke kinetic energy efficiency 
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inlet total-pressure recovery 
boundary -layer thickness 
boundary -layer displacement thickness 

Subscripts: 

C most forward point on cowl lip (see fig. 1(a)) 

CL 12° tangent point on cowl lip (see fig. 1(a)) 

t total 

th throat 

00 free stream 

1.0 mass -flow ratio of 1.0 (see appendix) 

Abbreviation: 

fms limiting mesh size 

ANALYSIS 
Analytical Method 

The analytical investigation that has been performed for the full-scale inlet of the 
HRE-AIM used an updated version of the computer program described in reference 5. 

The program may be exercised for either inviscid or combined viscid-inviscid solutions 
for real or perfect gas. Results obtained for this investigation were from real -gas, 
viscid-inviscid computations with laminar boundary -layer transition set to start at an X 
of 1.89 (corresponds closely to the transition regions of the HRE 2/3-scale inlet tests 
of ref. 1). 

The program incorporates several analytical methods that were listed in reference 5 
and are listed herein for convenience: 

(1) Method of characteristics 

(2) Blunt -body solutions (ref. 6) 


^R 

6 


6 


4 



(3) Sharp lip viscous interactions (ref, 7) 

(4) Laminar boundary layer (ref. 8) 

(5) Turbulent boundary layer (ref. 9) 

(6) Boundary -layer transition (momentum balance) 

(7) Shock— boundary -layer interaction (continuity) 

(8) Vortex sheets 

(9) Shock wave intersections 

Limitations of the computer program indicated in reference 5 are that it is valid for 
constant wall temperature only (requires some care in wall temperature selection to 
approximate heat transfer) and it does not compute shocks resulting from wave coales- 
cence. The blunt -body -method solution employed was noted to have convergence prob- 
lems in that trends for the pressure, density, velocity, and flow inclination yielded by the 
method are not as smooth as required for input to the method-of -characteristics solu- 
tions. However, the blunt -body -method solution is widely accepted and used in supersonic 
flow near blunt bodies. It was also noted in reference 5 that the computer solutions are 
sensitive to limiting mesh size (Zms) as are all method -of -characteristics computer 
programs. This and other sensitives encountered are discussed in the section "Results 
and Discussion." 


Analytical Model 

The model used in this analytical investigation was the full-scale inlet of the HRE- 
AIM which is shown schematically in figure 1(a). Coordinates of the inlet are listed in 
table 1, and the values for the cowl coordinates are for shock-on -lip position at Mach 6.0; 
other cowl lip locations were analyzed in which the appropriate cowl longitudinal coordi- 
nates were used. All X values are measured from the centerbody virtual vertex as 
shown in the insert sketch of the tip. The Xq values represent the distance from the 
vertex to the most forward point on the cowl leading edge as shown in the cowl insert. 
Note in the cowl detail the location of the most forward point and R(^) on the cowl 
lip in relation to the 12° tangent point ^Cl) used in reference 1. 

Throughout this investigation the X^ and R^ values are used; the Rq value is 
22.934 cm (9.029 in.). The inlet throat location indicated in figure 1(a) at X = 4.5 
(102.87 cm or 40.5 in.) was used in the present theoretical analysis. Surface meas- 
urements of the AIM have recently been reviewed, and the location of the throat station 
has been more accurately determined to be at 104.78 cm (41.25 in.). Change in the geo- 
metric location is not considered significant to the results because of the uncertainties in 
boundary -layer -thickness predictions and because the region of the throat displacement 
is one of nearly constant area. 
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Analytical Conditions 


Experimental tests of the HRE-AIM (photograph shown in fig. 1(b)) are scheduled 
to be conducted in the Lewis hypersonic tunnel facility (HTF); AIM is shown partially 
installed in the HTF in the photograph of figure 1(c). Analytical computations were there- 
fore performed in this investigation for the nominal conditions obtainable in the HTF which 
simulated high and low flight altitudes (see tables 2 and 3 and fig. 2). Calibration Mach 
numbers shown in table 2 were obtained during facility calibration tests discussed in ref- 
erence 10. The conditions actually used for the computer computations are listed in 
table 4. All 23 cases were computed by using the higher total pressures simulating the 
higher flight dynamic pressure conditions. However, for some Mach 5.0 cases at the 
higher total pressure, the characteristics solution terminated prior to the throat station; 
therefore, the cases at the lower total pressures were used as indicated in table 3 
(cases 1, 3, and 9). 


Method of Calculation 

Computer solutions are generated in two parts by the program. The external flow 
solution is first computed (generally referred to as job A) and the results in the region of 
the cowl lip are used to compute the internal flow (referred to as job B). These regions 
of flow are depicted in the sketch of figure 1(a). When a cowl lip shape and position are 
included in the input for running job A only, the program will also generate the cowl lip 
stagnation streamline from which the ratio of the captured free- stream mass flow and the 
additive drag coefficient are determined. 

RESULTS AND DISCUSSION 

The computer program cases listed in table 4 have been analyzed and are discussed 
in the following paragraphs. The results of these analyses are compared with some 
results which were presented in reference 1. 

Effect of Cowl Lip Positioning on Internal Flow 

Figures 3(a), 3(b), and 3(c) present typical internal flow shock patterns at various 
cowl lip positions for Mach numbers of 5.0, 6.0, and 7.0, respectively. The bottom sketch 
of each figure represents the case with the centerbody tip shock nearly impinging upon the 
cowl lip stagnation point. (For Mach 5.0 (fig. 3(a)) the cowl lip position in the bottom 
sketch is the same as for the Mach 6 shock-on-lip position.) From the bottom to the top 
of each page, sketches are arranged in order of increasing distance from the centerbody 
tip to the cowl lip. As this longitudinal distance increased, both the captured mass flow 
and the inlet throat area decreased. Also, the flow directed at the blunt cowl lip decreased 
in Mach number and increased in flow angle relative to the cowl lip internal surface 
because of greater external compression. This influenced the shape of the cowl lip shock 
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and the internal shock reflection patterns. These changes in the cowl lip shock shape and 
the internal shock reflection patterns, of course, affected the inlet performance, including 
pressure recovery and kinetic energy efficiency. It is interesting to note that for all Mach 
numbers and cowl positions considered, the cowl lip shock always impinged upon the center - 
body in the region of an X of 4.2. This similarity apparently results from a combination 
of the changes in the shape of the cowl lip shock and the inlet gap height at the cowl lip. 

A large region of wave coalescence occurred in the internal flow for all three Mach 
numbers and for the two most forward positions of the cowl. The computer program does 
not compute shocks that would result from wave coalescence such as shown in figure 3 
but instead employs a computing ’’bookkeeping" procedure generally accepted for the 
method of characteristics. The procedure is as follows: as a coalescence of two char- 
acteristic rays of the same family occurs, computation of the downstream ray is termi- 
nated and the remaining points of the previously computed upstream ray also become the 
remaining points of the downstream, ray. This procedure results in a discontinuity in the 
field of characteristics in the remaining computations, but without entropy and total - 
pressure losses that would be present if truly accounting for a new shock. For Mach 5.0 
and 6.0 (figs. 3(a) and 3(b), respectively), the large coalescence region exists only at the 
most forward cowl position (cases 2 and 11) but another region of coalescence occurs 
prior to the larger region for more rearward cowl positions. Then, as the rearward move- 
ment of the cowl continues, the large region is not apparent and only the small region 
emanating very close to the cowl lip remains. The large region of coalescence was also 
noted in reference 1 and was stated to be primarily the result of the compression fan 
originating from the centerbody aft compression surface (see fig. 1(a)) focusing on the 
cowl inner surface and reflecting from it. The small region of coalescence emanates 
from the cowl inner surface in the initial compression region along this surface. The 
absence of these small regions from the Mach 5.0 and 6.0 most forward positions of the 
cowl lip (cases 2 and 11) and the Mach 7.0 two most forward positions (cases 17 and 18) 
is not readily explainable. 

Schedules of Mass -Flow Ratio, Additive Drag Coefficient, 
and Aerodynamic Contraction Ratio 

Mass -flow ratio .- The computational results of the cases listed in table 4, along 
with the results of many job A computer cases, were used to establish the mass -flow- 
ratio curves presented in figure 4(a) which adequately define the ranges for the AIM tests. 
The solid lines represent mass flows for the three nominal test Mach numbers, 5, 6, 
and 7. Short -dash lines represent Mach numbers near the facility nozzle calibrated 
Mach numbers especially at Mach 5.0 and 7.0 (see table 2). The long-dash line at Mach 4 
is from reference 1. Values for the Mach 8 curve of reference 1 could not be obtained 
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for "spot check" computer cases; therefore, the Mach 8 curve of figure 4(a) was obtained 
from computations using job A. Values of the Mach 6 curve of figure 4(a) are in agree- 
ment with values presented in reference 1. Several cases were run with the centerbody 
tip shock near the cowl stagnation point to obtain accurate full -capture points (value of 
1.0) over the range from Mach 5.0 to Mach 8.0. At the bottom of figure 4(a) are small 
arrowheads depicting the presently planned cowl lip positions scheduled for testing during 
the AIM tests. 

Table 5 lists the points which define the curves for Mach numbers from 5.0 to 7.25 
of figure 4(a). These points were used to construct the curves presented in figure A1 of 
the appendix (the Mach 8 curve of fig. 4(a) was used to complete the figure up to Mach 8). 
Points from the curves of figures 4(a) and A1 were "curve fitted" and the mathematical 
expressions of the curve fits and their limitations are discussed in the appendix. Point 
values from the curves of figure 4(a) or some similar m;athematical expression is to be 
employed in an AIM data -reduction computer program. 

Additive drag coefficients .- Figure 4(b) contains the schedules of additive drag 
coefficients for Mach 5.0 to 7.25. The schedules correspond to the mass -flow -ratio 
schedules of figure 4(a). Values of additive drag coefficients which define the curves of 
figure 4(b) are given in table 6. These results may be used in "checking out" an AIM 
performance computer program. 

Contraction ratios .- The overall geometric -contraction -ratio schedule obtained 
from reference 1 is represented by the uppermost curve in figure 5. Using this schedule 
and the mass -flow -ratio values of figure 4(a), aerodynamic -contraction -ratio (geometric 
contraction ratio multiplied by the captured mass -flow ratio) schedules were obtained and 
are presented in figure 5. The aerodynamic -contraction -ratio curve for Mach 4 was 
obtained directly from reference 1. A faired curve through the peak values of the 
aerodynamic -contraction -ratio schedules is also depicted in figure 5; peak values are 
used in the discussion of the inlet total -pressure recovery. 

Throat Conditions 

Internal flow representation .- Results of the digital computer cases listed in table 4 
yielded the conditions for each case at the inlet throat station (X of 4.5). A typical inter- 
nal flow field produced from the computer results is shown in figure 6 with a tabulation of 
the conditions at the throat. The numbered points across the throat station are computed 
points in the characteristic net or are points obtained by straight-line interpolation 
between two net points. The throat profiles were determined from these point conditions. 

Throat profiles .- Throat profiles of local total-pressure ratio Pt/Pt oo and local 
Mach number are presented in figures 7 and 8. Figure 7 presents profiles for those cases 
where the spike tip bow shock is near the cowl lip stagnation point. (Note that at Mach 5.0 
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and 5.15 the cowl is located near shock-on-lip position at Mach 6.) Profiles for various 
cowl positions at Mach numbers of 5.0, 6.0, and 7.0 are presented in figures 8(a), 8(b), 
and 8(c), respectively. The discontinuities in the profiles are caused by the reflected 
internal shock crossing the throat station; a large discontinuity, however, was not always 
associated with the crossing of such a shock. Changes in the shapes of these profiles are 
attributed to the change of the internal shock patterns such as those shown in figure 3. 

The shapes of the total -pressure -ratio profiles appear to be similar to experimental 
results of the 2/3 -scale inlet investigation reported in reference 1, in that the total pres- 
sure is greatest near the centerbody. The lower total pressures next to the cowl are 
considered to be the result of the cowl lip bluntness affecting the flow next to the cowl 
surface. The spike tip bluntness effect is very insignificant in the flow next to the center - 
body at the throat because the effect on the flow is distributed about a larger area as the 
flow progresses downstream. 

Point values of the throat profiles in figures 7 and 8 represent the inviscid method- 
of -characteristics solutions. Similar throat condition profiles were generated by using 
characteristic mesh point conditions obtained from the reference 5 computer results as 
input to a computer program developed by one of the authors which determined mass - 
weighted average flow properties at the throat. The portions of the inviscid flow profiles 
within the boundary layers were replaced by using the boundary -layer -shape exponent 
values obtained from the reference 5 computer results and assuming a constant static 
pressure across the boundary layer equal to the pressure at the boundary -layer edge. 
(Superimposed in fig. 7 are values of the centerbody and cowl boundary -layer edge and 
displacement thicknesses for Mach 6.0 (case 11).) As a comparison, a total -pressure - 
ratio profile with the proper boundary -layer conditions is shown for Mach 6.0 full-scale 
theory from reference 1. Results of the mass -weighted averaging program included the 
total enthalpy losses and inlet performance parameters. 

Inlet Performance 

Total-pressure recovery .- Mass -weighted average total -pressure recoveries are 
presented in figure 9. Recoveries as a function of free -stream Mach number are shown 
in figure 9(a) and represent the results when the spike tip shock is nearly impinging on 
the cowl lip stagnation point. (Note that at Mach 5.0 the cowl is at the position for shock 
on lip at Mach 6.0.) The short -dash line in figure 9(a) represents theoretical mass- 
weighted results for the 2/3 -scale inlet T -model from reference 1 and the square symbols 
represent the corresponding mass-weighted 2/3-scale experimental data, also from ref- 
erence 1. The theoretical mass -weighted results from reference 1 for the full-scale 
inlet T -model are presented in figure 9(a) and do not exhibit good agreement with the 
2/3 -scale theoretical trend. (All the theoretical results were obtained by using the 
updated version of the ref. 5 computer program.) Direct comparison of the present AIM 
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theoretical results with the results from reference 1 (full-scale theory and the 2/3-scale 
theory and experimental data) can be made only for Mach 6.0; the agreement is best 
between the present AIM theory and the full-scale theory and is shown to be within 2 per- 
cent. The cowl positions of the present and the reference 1 Mach 6 full-scale cases are 
slightly different, which may contribute toward the difference of the recovery values. 

Pressure recoveries as a function of cowl position are presented in figure 9(b). 

The arrowheads represent the cowl positions planned for the AIM tests (also depicted in 
fig. 4(a)). Results for Mach 5.0 indicate a trend as represented by the solid line. The 
trend curve possesses slightly increasing pressure recovery values with increasing Xq 
values. The two largest values of cowl position have pressure recovery values that 
decrease slightly. These two cowl positions are on the decreasing portion of the 
aerodynamic -contraction -ratio curve (see fig. 5); therefore, the trend curve is shown to 
start decreasing at the Xq value corresponding to the peak aerodynamic contraction 
ratio. The trend of slightly increasing recovery values is also representative of two of 
the Mach 6.0 values ^at Xj;;; values of about 3.9 and 4.0?). At the point of the peak aero- 
dynamic contraction ratio for Mach 6.0 the decreasing trend is begun and it very nearly 
represents the last two cowl positions. A Mach 6.0 recovery value was obtained at an 
Xq of 4.0 which was low with respect to the trend. Computer results for this case 
indicated thick boundary layers internally (large ratios of boundary -layer thickness to 
throat channel height) that were about twice as thick as those for the other Mach 6.0 cases. 
Whether these thick boundary -layer results were truly valid could not be ascertained. 

The trend was not representative of the Mach 7.0 results; the results instead possessed 
a continuously decreasing trend. Recovery values for Mach 7.0, 7.15, and 7.25 at an X(^ 
of about 4.06 are for cases that possessed thick internal boundary layers which may 
account for these values appearing to be low with respect to the other Mach 7.0 recovery 
values. Since the trend curves are not representative of all cases, they should be treated 
with caution. 

Mach number .- Throat Mach numbers were obtained from the results of the mass- 
weighted averaging computer program previously mentioned. Throat Mach numbers as 
a function of free -stream Mach number are presented in figure 10(a). Points represent- 
ing the full-scale theory of reference 1 were obtained by averaging values of three char- 
acteristic net points near the throat station; this should result in slightly greater values 
than would be obtainable if enough points were given to account for the boundary -layer 
flow. The dash line is merely a fairing of the present AIM theoretical results. Note 
that the throat Mach number is 44 to 50 percent of the corresponding free -stream Mach 
number . 

Throat Mach numbers are shown as a function of the cowl lip position in figure 10(b). 
The straight-line trends are a family of constant slope (am/^Xq ~ -2.3^. At an Xq of 
3.747 for Mach 5.0, the value appears rather low because, as shown in the Mach number 


10 


profile (case 1) of figure 8(a), the Mach number drops off rapidly next to the centerbody, 
unlike the other profiles. The exact reason for this was not determined. The well- 
defined family of trends for the throat Mach number values are acceptable even though 
the pressure recovery values have poorly defined trends because the Mach number values 
are highly insensitive to total pressures and are very sensitive to static pressures. 

Kinetic energy efficiency .- Values of kinetic energy efficiency were obtained from 
the mass -weighted averaging computer program. The efficiencies are presented in fig- 
ure 11(a) as a function of free -stream Mach number for the cases with the spike tip shock 
nearly impinging on the cowl lip. (Note again that the Mach 5.0 and 5.15 cases have the 
cowl positioned for near shock on lip at Mach 6.0.) Efficiency values are shown in fig- 
ure 11(b) as a function of the cowl position. There is no readily evident trend of the 
efficiencies in either figure; the numerical average value for both figures is about 

0.93. 

The kinetic energy efficiencies were computed for a constant cooled-wall tempera- 
ture of 500 K (900° R) and resulted in throat total -enthalpy losses; ratios of the throat 
total enthalpy to the free -stream total enthalpy are listed in table 4. The reference 5 
computer program could not be operated for wall temperatures low enough to approximate 
the predicted wall temperatures of the AIM more closely. Enthalpy losses presently 
computed are therefore less than the expected losses for the AIM tests and thus the pres- 
ently computed kinetic energy efficiencies should be greater than the expected AIM test 
results. At Mach 6.0 the AIM inlet total -enthalpy losses are estimated to be about 9 per- 
cent, or a total -enthalpy -ratio value of 0.91. 

Effects of Computer Program Sensitivities 

Reference 5 indicated that the computer program employed for these analyses was 
sensitive to some of the input parameters. The effect of the mesh size and other sensi- 
tives noted during these analyses are discussed in the following paragraphs. 

Relationship between spike tip shock and cowl lip shock . - Computer difficulties were 
experienced when the spike tip shock was allowed to impinge near the cowl lip stagnation 
point. There was no reflection of the spike tip shock because the only possible solution 
in front of the blunt cowl lip is a normal shock wave. If the spike tip shock wave and the 
cowl lip shock wave were to combine in a reflection to the centerbody, a strong wave 
solution could exist. In one instance, for Mach 6.0, a resultant shock was computed that 
was a strong wave solution; however, a reflection of the shock from the centerbody never 
occurred. Attempts were also made to purposely allow the spike tip shock to impinge 
upon the cowl inner surface since such a test is planned for the AIM. These attempts 
were mostly unsuccessful except for one case where the spike tip shock did reach the 
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cowl inner surface after crossing the cowl lip shock but only a limited amount of flow- 
field computations was obtained; according to reference 5 this was, however, a complete 
solution. Successful complete solutions (internal shock reflections) were not obtained 
for the spike tip shock impinging near the cowl lip until the cowl location was changed 
slightly to allow the spike tip shock to impinge outside of the cowl lip stagnation point. 

(It should be noted in relation to the two foregoing shock problems that during the HRE 
2/3 -scale inlet model tests reported in ref. 1, the spike tip shock was focused as close 
as was physically possible to the cowl lip stagnation point with no detrimental efficiency 
losses. Also, the tip shock was allowed to enter the cowl without unstarting or drastically 
affecting the inlet operation.) 

Characteristic net limiting mesh size .- The limiting mesh size (Zms) has only slight 
effects on the accuracy of the results of the method of characteristics (the smaller, the 
more accurate); however, very often the selection of the correct Zms value was the 
determining factor of whether a complete and successful computation case was obtained. 
The standard ims value was 3.048 cm (1.2 in.), but in a few cases a more appropriate 
Zms value was required. 

A solution of interest was case 12 for Mach 6.0 which was run with three different 
limiting mesh sizes; the resulting internal flow shock patterns are presented in figure 12. 
With the largest limiting mesh size (4.572 cm (1.8 in.)) the first internal shock reflection 
occurred just downstream of the impinging shock. As the Zms value is decreased, the 
reflection occurs farther downstream. A limiting mesh size of 1.524 cm (0.6 in.) 
resulted in the shock reflecting downstream of the throat section. The trend of shifting 
the shock reflection downstream with decreasing limiting mesh size appears to be the 
result of the cowl lip shock intersecting the centerbody boundary layer in an increasingly 
more asymptotic manner; the trend is not definite, however, as it was for case 18 
(fig. 3(c)) at Mach 7.0. When limiting mesh size was decreased for this case, the reflec- 
tion of the shock was shifted upstream. It appears that the mesh size changes the shape 
of the cowl lip shock in that the mesh size affects the boundary-layer growth along the 
centerbody and thus affects the flow conditions approaching the cowl lip. However, no 
definite trend of the effect of the limiting mesh size upon the solutions can be concluded. 
Along the same lines, the boundary -layer solution is sensitive to the surface contour and 
to the resulting relative locations of expansions and shock wave reflections which, as 
discussed, are affected by the limiting mesh size. For cases 17 and 18 at Mach 7.0 
(fig. 3(c)), the edge of the boundary layer and the displacement thickness appeared to be 
distorted. This distortion is thought to be the result of the relative location of the expan- 
sions and reflected shock waves. A more appropriate Zms value may resolve the dis- 
tortions; however, several different values were attempted. 
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Even though this computer program has limitations and parts which could be 
improved, the program as a whole appears to be the best available for such inlet design 
analyses as are presented here. 


SUMMARY OF RESULTS 

A theoretical full-scale inlet investigation has been conducted in support of the 
technology development of the NASA Hypersonic Research Engine (HRE) Project which 
is scheduled to be culminated with the experimental investigation using the Aerothermo- 
dynamic Integration Model (AIM). The theoretical analysis used a method -of - 
characteristics computer program for defining the AIM inlet flow conditions (real -gas, 
viscid-inviscid flow) at various inlet cowl lip positions over the range of AIM free- 
stream test conditions at Mach numbers from 5.0 to 7.25. 

Mass -flow ratios and additive drag coefficients were well defined over the range of 
AIM test conditions. The mass -flow values and mathematical curve fitting expressions 
obtained for these values are suitable for employment in the AIM data -reduction program. 

Mass -weighted average inlet performance parameters, including total -pressure 
recovery, throat Mach number, and kinetic energy efficiency, as functions of free -stream 
Mach number and cowl positions were obtained. The results for these performance 
parameters can be summarized as follows: 

1. The pressure recovery decreases with free -stream Mach number. 

2. At a free -stream Mach number of 5.0 the pressure recovery increased linearly 
for a range of cowl positions. 

3. The throat Mach number decreased linearly as a function of cowl position for 
various free -stream Mach numbers. 

4. Kinetic energy efficiency results as functions of both free -stream Mach number 
and cowl position did not possess any definite trends; however, a numerical average of 
all kinetic -energy -efficiency values was about 0.93. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., Jime 11, 1974. 
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APPENDIX 


CURVE -FIT EXPRESSIONS 

Cowl Positions as a Function of Mach Number for 
Constant Values of Mass-Flow Ratios 

Quadratic curve fits for points from the family of curves presented in figure A1 have 
been obtained by using a small digital computer program. Two "fits" were obtained for 
each curve; one for the Mach number limit from 5.0 to 6.0 and the other for Mach num- 
bers from 6.0 to 7.25. Table A1 presents computer point value results for each mass- 
flow-ratio curve. The quadratic expression of the curve fits are shown at the bottom of 
table Al, just below the tabulation of the proper values for constants A, B, C, D, E, and F. 
These point values are in reasonably good agreement with the curves of figure Al. 



Cowl position, 
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Figure Al.- Effect of free-stream Mach number and cowl position on constant values 

of mass-flow ratios. 



TABLE Al.- COMPUTER POINT VALUE RESULTS FROM QUADRATIC EXPRESSIONS OF CURVES IN FIGURE A1 
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APPENDIX 


Mass- Flow Ratios as Functions of Mach Number and Cowl Position 


An exponential curve fit was obtained of points defining the mass -flow ratios; the 
mass-flow ratios were approximated reasonably well with certain limitations. The 
expression for the curve fit was obtained by using the difference between the curve for 
a mass-flow ratio of 1.0 and the other curves of figure Al. The difference 1.0 - 

Xc - Xc 1 0 

was plotted in figure A2 as a function of a cowl position parameter — ^ for the 


xc,i.o 


/V1oq = 8. 0-] 


8 



AX = 


1.0 


''C. 1.0 

Figure A2.- Mass-flow parameter. 

various free-stream Mach numbers. The values for each Mach number are connected 
with a straight-line fairing which reasonably well represents the point values in the 

range of 1.0 - from 0.1 to 0.3. The expression for such curves is 

^ 00 
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APPENDIX 


y = ae2*3026bx 

where, for the present curves, 


(Al) 


y = 1-0 - 

^ nioo 


Xp - Xp 1 Q 
X = Ax = 5=^^ 


Xc,1.0 

a represents the y -intercept, and b is the slope expressed as 




(y^c - xc,i,o ^ 

V ^c,i.o , 


For 5.0 £ M ^ 6.0 

Xc 1 0 = -0.213935 + 1.130883M - 0.074686 m2 

and for 6.0 £ M £ 7.25 

Xc 1 0 = 0.665611 + 0.866032M - 0.054976 m2 

A plot of the y -intercepts as a function of free -stream Mach number is somewhat 
random and is therefore represented by a straight line in figure A3; the equation for this 
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Figure A3. - y -intercepts of the mass -flow parameters. 
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APPENDIX 


straight line defines the y -intercepts for the various Mach numbers as 
a = 0.001333M + 0.019334 

Values of the slopes determined from figure A2 were plotted in figure A4 as a function 
of Mach number. 



Figure A4,- Slopes of the mass -flow parameters. 

A quadratic expression was obtained for the curve from Mach 5.0 to 6.0 and another 
was obtained for the curve from Mach 6.0 to 7.25. The expression for 5.0 = M = 6.0 is 

b = Bj = 126.30424 - 51.12572M + 5.51528M^ 

and for 6.0 = M = 7.25 is 

b = B2 = 105.899902 - 42.936332M + 4.717169 m2 
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APPENDIX 


Therefore, by substituting the expressions just discussed into equation (Al), the 
mass-flow ratio is defined as 


^ = 1.0 - ae2-3026bx 
nioo 


(A2) 


Equation (A2) was used in developing a small digital computer program which produced 
results that agreed with the curves of figure A2, that is, matched the point values as 
closely as the straight line matched the values. The mass-flow ratio was forced to a 
limiting value of 1.0. 
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TABLE 1.- HRE-AIM INLET COORDINATES USED IN THE 


PRESENT AIM THEORETICAL ANALYSIS 
jMach 6.0 cowl design position] 


(a) Spike 


(b) Cowl 


X = 0.066 
R= 0 



X 

1 

R 

0.066 

0.0 

.077 

.014^ Straight Line 

2.040 

o 

o 

o 

CD 

CO 

2.145 

.379 

2.271 

.404 

2.410 

.432 

2.537 

.458 

2.650 

.482 

2.875 

.531 

2.974 

.554 

3.100 

.584 

3.212 

.613 

3.295 

. 6 36^^ Straight Line 

3.373 

.658J 15.819° 

3.640 

.740 

3.787 

.793"! Straight Line 

4.190 

.956j 22.0° 

4.230 

.971 

4.274 

.984 

4.308 

.991 

4.341 

.997 

4.374 

1.003 

4.408 

1.087 

4.411 

1.013 

4.500 

1.020 Throat 

4.602 

1.030 

4.659 

1.035 

4.714 

1.040 



X 

R 

3.876 

1.0 

12° 

3.933 

1.012 


3.986 

1.021 

lOO 

4.019 

1.027 


4.046 

1.031 

8° 

4.085 

1.036 


4.166 

1.044'| Straight Line 

5.326 

1.159J 

5.645° 




TABLE 2,- EXPERIMENTAL FREE-STREAM CONDITIONS 


Nominal 

Mach 

number 

Pt, 

oc 

Tt 

,00 

Simulated altitude 
(a) 

Calibrated 

Mach 

number 

(b) 

atm 

psia 

K 

Or 

km 

ft 

5.0 

14.08 

207 

1248 

2245 

25.30 

83 000 

5.18 

5.0 

28.57 

420 

1225 

2204 

20.71 

68 000 

5.18 

6,0 

31.70 

466 

1665 

2995 

27.72 

91 000 

6.05 

6,0 

63.22 

930 

1628 

2930 

23.16 

76 000 

6.05 

7.0 

^^34. 00 

^500 

2211 

3980 

34.59 

113 500 

7.25 

7.0 

^68.00 

^1000 

2130 

3835 

29.87 

98 000 

7.25 


^For nominal Mach numbers and free-stream static pressures, 
*^From reference 10. 

^Arbitrarily selected as half of AIM design pressure. 
^Maximum AIM design operating pressure. 


TABLE 3.- HRE FLIGHT ENVIRONMENT 


Mach 

number 

Altitude 

Static 

pressure, 

Poo 

Static 

temperature, 

T 

C30 

Total 

pressure, 

Pt,^ 

Total 

temperature, 

'Pt,OC 

Dynamic 

pressure, 

^OC 


km 

ft 

atm 

psia 

K 

OR 

atm 

psia 

K 

OR 

atm 

psfa 

5.0 

20.11 

66 000 

0.0536 

0.788 

216.8 

390.1 

30.95 

455 

1222 

2200 

0.937 

1983 


20.71 

68 000 

.0481 

.707 

217.3 

391.2 

28.23 

415 

1225 

2204 

.852 

1802 


25.3 

83 000 

.02401 

.353 

222.0 

399.3 

14.08 

207 

1248 

2245 

.420 

889 

6.0 

22.24 

73 000 

.03842 

.565 

218.9 

393.9 

72.75 

1070 

1628 

2930 

.964 

2041 


23.16 

76 000 

.03333 

.490 

219.8 

395.5 

63.22 

930 

1628 

2930 

.840 

1777 


27.72 

91 000 

.0166 

.244 

224.2 

403.7 

31.70 

466 

1665 

2995 

.4175 

884 

7.0 

24.38 

80 000 

.02769 

.407 

221.0 

397.7 

154.5 

2270 

2089 

3760 

.9475 

2005 


25.17 

82 600 

.0245 

.360 

221.8 

399.1 

137.5 

2020 
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3770 

.840 

1776 


29.72 

97 600 

.01259 

.185 

225.0 

406.9 

70.0 

1030 

2129 

3830 

.4322 

915 
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TABLE 4,- CONDITIONS FOR AND RESULTS FROM THE THEORETICAL COMPUTER PROGRAM CASES 
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incomplete throat flow conditions. 



TABLE 5.- VALUES DEFINING THE CURVES^ OF MASS-FLOW RATIOS AS A 
FUNCTION OF COWL POSITION FOR CONSTANT MACH NUMBERS 



^See figures 4(a) and Al. 
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TABLE 6.- VALUES DEFINING THE CURVES^ OF ADDITIVE DRAG 
COEFFICIENTS AS A FUNCTION OF COWL POSITION FOR 
CONSTANT MACH NUMBERS 


Cowl lip position Xq at M 


3.5733 

3.7780 

3.8610 

3.9120 

3.9530 

3.9835 

4.0095 

4.0300 

4.0450 

4.0575 

4.0660 

4.0730 

4.0815 

4.0900 

4.0970 

4.1035 

4.1100 

4.1185 

4.1255 

4.1325 

4.1395 


3.6227 

3.8070 

3.8860 

3.9340 

3,9720 

4.0010 

4.0240 

4.0415 

4.0555 

4.0645 

4.0725 

4.0800 

4.0875 

4.0950 

4.1015 

4.1090 

4.1166 

4.1230 

4.1300 

4.1370 

4.1425 


3.8827 

3.9535 

3.9975 

4.0280 

4.0505 

4.0670 

4.0800 

4.0895 

4.0965 

4.1030 

4.1110 

4.1180 

4.1240 

4.1300 

4.1350 

4.1400 

4.1445 

4.1500 



4.0830 
4.0925 
4.1000 
1080 
1145 
1210 



1430 
4.1480 
4.1525 
4.1575 
4.1620 
4.1655 


4.0307 

4.0615 

4.0820 

4.0965 

4.1075 

4.1150 

4.1200 

4.1235 

4.1280 

4.1315 

4.1350 

4.1390 

4.1425 

4.1465 

4.1500 

4.1540 

4.1580 

4.1615 

4.1650 

4.1690 

4.1730 


4.0453 
4.0695 
4.0880 
4.1005 
4.1115 
4.1185 
4.1230 
4.1275 
4.1305 
4.1335 
4.1380 
4.1415 
4.1455 
4.1495 
4.1530 
4.1570 
4.1505 
4.1650 
4.1680 
4.1720 
4.1755 



4.1475 

4.1510 

4.1550 

4.1590 

4.1525 

4.1665 

4.1695 

4.1735 

4.1780 


See figure 4(b). 
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(a) AIM schematic. 

Figure 1.- Hypersonic Research Engine (HRE) Aerothermodynamic Integration Model (AIM) 




(b) Pretest photograph of AIM 
Figure 1.- Continued. 
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(c) AIM partially installed in Lewis hypersonic tunnel facility (HTF). 

Figure 1.- Concluded. 



Static-pressure altitude, 



5 6 7 8 


Free-stream Mach number, 

Figure 2.- Free-stream conditions. (Altitude values obtained by using static - 
pressure values corresponding to total pressures and Mach numbers in 
tables 2 and 3.) 
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(a) Mach 5.0. 

Figure 3.- Internal flow -field shock patterns for various cowl positions 


























































Figure 8.- Continued. 
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Figure 11.- Inlet mass -weighted average kinetic energy efficiencies, present AIM theory. 
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Figure 12.- Ef: 



